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We address the mechanism of drag reduction by polymers in turbulent wall bounded flows. On the 
basis of the equations of fluid mechanics we present a quantitative derivation of the "maximum drag 
reduction (MDR) asymptote" which is the maximum drag reduction attained by polymers. Based 
on Newtonian information only we prove the existence of drag reduction, and with one experimental 
parameter we reach a quantitative agreement with the experimental measurements. 



The addition of few tens of parts per million (by 
weight) of long-chain polymers to turbulent fluid flows 
in channels or pipes can bring about a reduction of the 
friction drag by up to 80% 0, U H • This phenomenon 
of "drag reduction" is well documented and is used in 
technological applications from fire engines (allowing a 
water jet to reach high floors) to oil pipes. In spite of a 
large amount of experimental and simulational data, the 
fundamental mechanism has remained under debate for a 
long time 0,0,0- I n such wall-bounded turbulence, the 
drag is caused by momentum dissipation at the walls. For 
Newtonian flows (in which the kinematic viscosity is con- 
stant) the momentum flux is dominated by the so-called 
Reynolds stress, leading to a logarithmic (von-Karman) 
dependence of the mean velocity on the distance from 
the wall 0]. However, with polymers, the drag reduction 
entails a change in the von-Karman log law such that a 
much higher mean velocity is achieved. In particular, for 
high concentrations of polymers, a regime of maximum 
drag reduction is attained (the "MDR asymptote"), in- 
dependent of the chemical identity of the polymer 0, 
see Fig. 1. In this Letter we elucidate the fundamen- 
tal mechanism for this phenomenon: while momentum is 
produced at a fixed rate by the forcing, polymer stretch- 
ing results in a suppression of the momentum flux from 
the bulk to the wall. Accordingly the mean velocity in 
the channel must increase. We derive a new logarithmic 
law for the mean velocity with a slope that fits existing 
numerical and experimental data. The law is universal, 
thus explaining the MDR asymptote. 

Turbulent flows in a channel are conveniently discussed 
for fixed pressure gradients p' = —dp/dx where x, y 
and z are the lengthwise, wall-normal and spanwise di- 
rections respectively. The length and width of the chan- 
nel are usually taken much larger than the mid-channel 
height L, making the latter a natural re-scaling length 
for the introduction of dimensionless (similarity) vari- 
ables. Thus the Reynolds number Tie, the normalized 
distance from the wall y + and the normalized mean ve- 
locity V + (y + ) (which is in the x direction with a depen- 
dence on y only) are defined by 

Tie = L^p7L/u , y+ = yRejL , V+ = V/^L , (1) 

where vq is the kinematic viscosity. One of the most 
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FIG. 1: Mean normalized velocity profiles as a function of 
the normalized distance from the wall during drag reduction. 
The data points from numerical simulations (green circles) 
[lH and the experimental points (open circles) |l2jj represent 
the Newtonian results. The black solid line is our theory 
Eg. 1131 . which for large y + agrees with von-Karman's loga- 
rithmic law of the wall (|5J . The red data points (squares) [l3i| 
represent the Maximum Drag Reduction (MDR) asymptote. 
The dashed red curve represents our theory for the profile 11911 
which for large y + agrees with the universal law l|16jl . The 
arrow marks the crossover from the viscous linear law 1121 to 
the asymptotic logarithmic la w lllrit . The blue filled triangles 
[l3^ and green open triangles |14| represent the cross over, for 
intermediate concentrations of the polymer, from the MDR 
asymptote to the Newtonian plug. Our theory is not detailed 
enough to capture this cross over properly. 



famous universal aspects of Newtonian turbulent channel 
flows is the "logarithmic law of the wall" which in these 
coordinates is expressed as 



V+(y + ) 
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for y + > 30 , (2) 



where the von-Karman constant n K ss 0.436 and inter- 
cept B re 6.13 [Ig . For y+ < 10 one observes a viscous 
sub-layer, V + (y + ) = y + , Fig. Q The riddle of drag re- 
duction is then introduced in relation to this universal 
law: in the presence of long chain polymers the mean ve- 
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locity profile V + (y + ) (for a fixed value of p' and channel 
geometry) changes dramatically. For sufficiently large 
concentration of polymers V + (y + ) saturates to a new 
(universal, polymer independent) "law of the wall" 0, 



V + (y + ) = 11.7 lny+ - 17, for y+ > 12 



(3) 



For smaller concentration of polymers the situation is as 
shown in Fig. The Newtonian law of the wall (J2J) is the 
black solid line for y+ > 30. The MDR asymptote © is 
the dashed red line. For intermediate concentrations the 
mean velocity profile starts along the asymptotic law JSJ, 
and then crosses over to the so called "Newtonian plug" 
with a Newtonian logarithmic slope identical to the in- 
verse of von-Karman's constant. The region of values of 
y + in which the asymptotic law @ prevails was termed 
"the elastic sub-layer" 0. The relative increase of the 
mean velocity (for a given p') due to the existence of 
the new law of the wall J2J is the phenomenon of drag 
reduction. Thus the main theoretical challenge is to un- 
derstand the origin of the new law ©, and in particular 
its universality, or independence of the polymer used. A 
secondary challenge is to understand the concentration 
dependent cross over back to the Newtonian plug. In this 
Letter we argue that the phenomenon can be understood 
mainly by the influence of the polymer stretching on the 
y + -dependent effective viscosity. The latter becomes a 
crucial agent in carrying the momentum flux from the 
bulk of the channel to the walls (where the momentum is 
dissipated by friction). In the Newtonian case the viscos- 
ity has a negligible role in carrying the momentum flux; 
this difference gives rise to the change of Eq. (J2J) in favor 
of Eq. which we derive below. 

The equations of motion of polymer solutions can be 
written as 



(4) 



dU/dt + U ■ VU = -Vp + V • T + is V 2 U 



where T is the extra stress tensor that is due to the poly- 
mer. Denoting the polymer end-to-end vector distance 
(normalized by its equilibrium value) as r, the average 
dimensionless extension tensor 7Z- is IZij = (r^j), and 
the extra stress tensor is (with = dUi/dxj), 

T = u p {uj ■ Tl + 7L ■ uj t - 871/ dt - U ■ VR.) . (5) 

Here v v (proportional to the polymer concentration) is 
the polymeric contribution to the viscosity in the limit 
of zero shear. In order to develop a transparent the- 
ory we propose to ignore the fluctuations of 7Z- compared 
to its mean. In other words, we will take 1Z- ~ (7V). 
Simplifying further the tensor structure, assuming that 
IZij = IZSij we can rewrite Eq. Q| in the form 

dUi/dt + UjVjUi = -ViP + VjV (wy + Uji) , (6) 
v = v + v p n, P = p + du/dt + U -Vv , (7) 

with a scalar (y-dependent) effective viscosity v. 



Armed with the effective equation we establish the 
mechanism of drag reduction following the standard 
strategy of Reynolds, considering the fluid velocity U(r) 
as a sum of its average (over time) and a fluctuating part: 

U(r,t) = V(y)+u(r,t) , V(y) = (U(r,t)) . (8) 

For a channel of large length and width all the averages, 
and in particular V(y) =>■ V(y), are functions of y only. 
The objects that enter the theory are the mean shear 
S(y), the Reynolds stress W(y) and the kinetic energy 
K(y); these are defined respectively as 

S(y) = dV{y)/dy , W(y) = -(u x u y ) , K(y) = (\u\ 2 )/2. 

A well known exact relation y| between these objects is 
the point-wise balance equation for the flux of mechanical 
momentum; near the wall (for y <C L) it reads : 

v(y)S(y) + W(y) = p'L . (9) 

On the RHS of this equation we see the production of mo- 
mentum flux due to the pressure gradient; on the LHS 
we have the Reynolds stress and the viscous contribu- 
tion to the momentum flux, with the latter being usually 
negligible (in Newtonian turbulence v = v§) everywhere 
except in the viscous boundary layer. The y dependence 
of the effective viscosity v(y) in the elastic layer will be 
shown to be crucial for drag reduction. 

A second relation between S(y), W(y) and K(y) is 
obtained from the energy balance. The energy is created 
by the large scale motions at a rate of W(y)S(y). It is 
cascaded down the scales by a flux of energy, and is finally 
dissipated at a rate e, where e = iy(uj^). We cannot 
calculate e exactly, but we can estimate it rather well at 
a point y away from the wall. When viscous effects are 
dominant, this term is estimated as v{a/y) 2 K(y) (the 
velocity is then rather smooth, the gradient exists and 
can be estimated by the typical velocity at y over the 
distance from the wall). Here a is a constant of the order 
of unity. When the Reynolds number is large, the viscous 
dissipation is the same as the turbulent energy flux down 
the scales, which can be estimated as K(y)/T(y) where 
r(y) is the typical eddy turn over time at y. The latter 
is estimated as y/b^jK(y) where b is another constant 
of the order of unity. We can thus write the balance 
equation at point y as 

V {y){a/yfK{y) + bK^ 2 (y)/y = W(y)S{y) , (10) 

where the bigger of the two terms on the LHS should win. 
We note that contrary to Eq. © which is exact, Ea. H10|) 
is not exact. We expect it however to yield good order of 
magnitude estimates as is demonstrated below. Finally, 
we quote the experimental fact 0, 0] that outside the 
viscous boundary layer 



W(y) 
K(y) 



for Newtonian flow, 



for viscoelastic flow. 



(11) 
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The coefficients c N and c v are bounded from above by 
unity. (The proof is \c\ = \W\/K < 2\(u x u y }\/ (u 2 x + 
Uy) < 1, because (u x ± u y ) 2 > 0). While c N is known 
quite accurately, c N 0.5, the actual values of c v varies 
somewhat from experiment to experiment, such that the 
ratio c v /c N (which is all that we need to use below) varies 
between 0.3 and 0.7. In our estimates we take c v /c N = 
0.5. 

We show now that Eqs. (0 - lllfl are sufficient for 
deriving the Newtonian law © and the viscoelastic law 
with equal ease. Begin with the Newtonian case. The 
result of this derivation is not new - but we want to stress 
that the same equations give rise to both the well known 
and the new results. Substitute Eq. l(TT|) in Eqs. (P HU|) 
with v(y) = turning them into algebraic equations for 
K(y) and S(y), and eventually to a first order differential 
equation for V(y). In the viscous sub-layer K{y) = 0, 
and the solution in re-scaled coordinates is 

V+(y + )=y + , y + <y+ = a/c n , (12) 

Outside the viscous sub-layer (y + > ?/+) we find 

V + (y + ) = k- 1 \nY(y+) +B- A(y+) , (13a) 
B = 2y+ - n- 1 In [e (1 + 2 k k j/+)/4k k ] , (13b) 



where k k = c N /6 and we defined 



Y{y + ) 
A(y + ) 



+ 2 + 
/+ ~Vv 



'2, 



2nlyf +Ak k [Y{ V +)- V +] +1 
2nly+ 



For y + > y+ Eqs. {T3J) turns into © since Y(y + ) — ► y+ 
and A(y+) -► 0. 

Note that Eqs (|13fl pertain to the whole y + domain. By 
taking the experimental values of k k and B we compute 
yy ~ 5.6 to be compared with the experimental value of 
5.5±0.5, cf. Q. The resulting Newtonian profile, Eqs. 
(|12fl and (T3J, is shown in Fig. 1 as the black solid line. 
The excellent agreement with the experimental and nu- 
merical data in the entire region of y + indicates that our 
balance equations are sufficiently accurate and we can 
proceed to the viscoelastic case. 

To see how the law © emerges we consider Eqs. © 
and i|10|) with y-dependent effective viscosity. We should 
warn the reader that this is not fully justified - there 
are terms in the full viscoelastic equations which cannot 
be simplified to the form of a space dependent viscos- 
ity. Nevertheless we propose (and justify further below) 
that the terms with effective viscosity are the main terms 
that allow the momentum flux to be carried in the elas- 
tic sublayer. In other words, when the concentration 
of the polymer is large enough we can neglect the sec- 
ond term in favor of the first in Eq. @ and estimate 
v{y) = Lp'/S(y). Substituting this estimate in Eq. (|10f> . 



neglecting the second term on the LHS, using Eq. Ijllll . 
and finally re-scaling the variables results in 



dV + /dy + = c N y+/c v y H 



(14) 



From here follows immediately the new "logarithmic law 
of the wall" 



V+(y + ) = K v - 1 lny+ + S v 



C v/ C N^' ( 15 ) 



where the intercept B v is still unknown (but will be de- 
termined momentarily). The slope of the new law is inde- 
pendent of the polymer concentration and is greater than 
the von-Karman slope 1/k k ; the ratio of the slopes is 
in fact K K y+c N /c v , which is about 4.9 according to the 
above estimates c N /c v w 0.5 and y + w 5.6. Compar- 
ing with the measured ratio of slopes in Eqs. (J2J and 
© which is about 5.1, we consider our estimates to be 
quite on the mark. As said before, this increase in slope 
is the phenomenon of drag reduction. We stress that 
the information gained from the Newtonian data alone is 
sufficient to predict drag reduction, since c v < 1. 

To find B v we match the logarithmic law i|15|) to the 
viscous sublayer solution (|12[) by the value of V + (y + ) and 
its derivative. First, the logarithmic law has slope 1 at 
y+ = n~ l . Next, matching at this point the viscous so- 
lution V + {y^ n ) = k^ 1 to fTB^I we find B v — ln(e k v )/k v . 
We note that if we substitute the experimental value 
k^ 1 = 11.7 we find B v = — 17 in perfect agreement with 
Eq. J2J. We thus write the law Q15(l in its final form 
(with just one constant remaining) 



V+(y+) = —\n(eK v y+) . 



(16) 



Note that this universal result is obtained without ref- 
erence to any model of the polymer dynamics, and the 
only assumptions are that the polymer viscosity domi- 
nates the momentum transfer in the elastic layer and the 
logarithmic law (|15J) is valid all the way to y+. At this 
point we demonstrate that these two assumptions, which 
lead to the universality of ifTfijl . are well supported by a 
closer consideration of the polymer dynamics. 

The fundamental reason for the appearance of a y- 
dependent viscosity is the coil-stretch transition of the 
polymers under turbulent shear. This transition occurs 
in any reasonable model of the polymer-turbulence inter- 
action, the simplest of which is of the form 



dTl/dt = -K/T p + s(y)K , s(y) = y^4), (17) 

where f p is the polymer relaxation time. In the elastic 
sublayer we estimate s(y) — gy/K{y)/y, with some con- 
stant g w O(l), and expect a coil stretch transition when 
g\/K(y)/y = l/fp. Defining r p = gf p we then find in 
the elastic sublayer (where polymers are stretched) 



K{y) = (y/r p y 



(18) 
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This equation is important since it can be solved together 
with Eqs. © and IjlOfl to find S + (y + ) when the polymers 
are stretched. The answer is 

+ ^ + [(2» K ,^,) 2 -l](^) 2 }. 



Navier-Stokes equations with a space dependent viscos- 
ity according to the above theory. We predict that a 
viscosity profile that remains constant for < y + < y+, 
and then grows linearly towards the center, should result 
in drag reduction in much the same way as seen in ex- 
periments with polymers. Such simulations would render 
direct support to the views offered in this Letter. 



where £ — 2 K K y+ c v r p y/ hp' . Equation H14() is recaptured 
from this, more general equation, when Tie is large and 
{...}—► 1. Then the identity of the polymer is lost, giv- 
ing rise to the universal drag reduction asymptote 
It can be also checked that the matching point y+ used 
above indeed connects smoothly the viscous sublayer to 
the asymptotic law (|16J) (see arrow in Fig. 1). We can 
therefore conclude that our derivation of Eq. 1(16(1 is fully 
consistent with the polymer dynamics, and we under- 
stand how the polymer characteristics drop out, leading 
to the universal law (|16ll . 

The mechanism of drag reduction is then the suppres- 
sion of the Reynolds stress in the elastic sublayer. The 
Reynolds stress is the main agent for momentum flux in 
the Newtonian case, and its suppression results in an in- 
crease of the mean mechanical momentum (velocity) in 
the channel. The increase in viscosity of course leads to 
increased energy dissipation, but this is immaterial for 
the phenomenon. This is the main difference between 
drag reduction in wall-bounded turbulence and in homo- 
geneous turbulence [3] ■ We note that the theory predicts 
a very low value of K(y) in the elastic layer. In reality in 
every wall bounded turbulent flows there will be ejecta of 
parcels of fluid with high level of K(y) from the Newto- 
nian region towards the walls. These would be picked up 
in measurements, and would seem to contradict the low 
value of K (y) predicted above. In fact there is no real 
contradiction since such random ejecta do not contribute 
much to the momentum flux that is so crucial to our dis- 
cussion. We also note a term in the dynamical equations 
for the viscoelastic flow of the form d(rirf) jdt, which can- 
not be written as an effective viscosity. This term van- 
ishes for uncorrelated rotations of individual polymeric 
molecules. However it can contribute to the energy flux 
from the mean shear to turbulent fluctuations in wall 
bounded flows due to a possible "correlation instability" , 
which leads to synchronized rotation of neighboring poly- 
mers. A more detailed consideration of such terms should 
provide also the cross over behavior seen in Fig. 1. 

Finally we propose direct numerical simulations in 
channel flows to test our approach. Instead of simulat- 
ing the viscoelastic equations, we propose to simulate the 
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